RADIO SCIENCE Journal of Research NBS/USNC-URSI 
Vol. 68D, No. 2, February 1964 



Diurnal Changes of Transmission Time in the Arctic 
Propagation of VLF Waves 



W. T. Blackband 

Contribution from the Royal Aircraft Establishment, Farnborough, Hants, England 

(Received September 10, 1963) 

The diurnal variation of the transmission time for the signal from a distant VLF trans- 
mitter results from the solar illumination of the lower ionosphere. For a path crossing the 
Arctic regions there will be no diurnal change for that part of the path which is illuminated 
by the midnight sun. It is shown that for a path crossing the Arctic Circle the diurnal 
change curve is of the normal trapesium shape at the equinoxes but that it takes on different 
forms at midsummer and midwinter. An analysis of experimental curves shows that they 
are of the form predicted. A simple rule for computing the change over on Arctic path is 
shown to agree well with the experimental data available. 



1. Introduction 

The diurnal variations in the transmission time 
for signals received at VLF have been studied in 
some detail for transmission paths in the lower 
latitudes but propagation in the Arctic regions has 
not been studied as closely. The two factors which 
introduce complexity into the changes over Arctic 
paths are 

(a) the periods of 24 hr daylight or darkness at 

certain seasons, 

(b) the large differences in local time between the 

ends of relatively short great circle paths. 

During the midsummer period of 1961 experi- 
mental observations of the diurnal variations in 
transmission time have been made at two points 
within the Arctic Circle, Bod0 in Norway and 
Sondrestrom in Greenland, and also in Farnborough 
for the path to Seattle which lies partly within the 
Arctic Circle. 

In this paper an explanation is given of the 
mechanism of these changes, and predicted changes 
are compared with those observed. The experi- 
mental data relate to frequencies in the VLF band 
between 16 and 20 kc/s. 

2. Diurnal Change in Phase at Middle 
Latitudes 

This change is brought about by the change in 
the apparent height of the ionosphere. For rays 
near normal incidence this change is from about 
90 kiii. in darkness to about 72 km in full daylight, 
near glancing incidence these heights are somewhat 
less. In figure 1 is shown a ray from a transmitter 
T reflected from the ionosphere to a receiver R. If 
the receiver is west of the transmitter, then when 
it is midnight at T, the path will all be in darkness 
and the ionosphere will be at its full height of 90 km. 
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Figure 1. Reflections of signal from the ionosphere by day and 
by night. 



Between the times of dawn at the receiver (D R ) and 
sunset at the transmitter (S T ), the whole path is in 
daylight and the ionosphere is at its lower height 
of about 72 km. 

During the daytime the ray between T and R 
will follow the dotted path of figure 1 which is shorter 
than the path up to the nighttime reflecting level 
which is shown with a full line. This reduction in 
path length reduces the transit time of the ray by 
At, this change being apparent as an advance in the 
phase of the received signal. While the value of At 
over a given path is largely independent of frequency, 
the phase change varies almost linearly with frequency 
and for this reason the diurnal phase changes observed 
at different frequencies are most conveniently com- 
pared on the basis of their equivalent time delays. 

Between the periods of all darkness and all light 
along the path TR there are two periods of changing 
state when the transit time has intermediate values. 
These are shown by the sloping parts of the graph 
of figure 2. At ranges above 4500 km where mode 
interference cannot complicate the pattern, the 
joins of the two steady states are approximately 
straight with fairly definite angles to the "diurnal 
change trapezium" — D T l) R S T S R . 

The variation of At with length of transmission 
path lias been discussed previously [Blackband, 
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Figure 2. The diurnal change t? % apezium. 



1961 and 1963]. A curve based on measured data 
is drawn with a solid line in figure 3. 1 In this figure a 
dotted line shows the variation of At against distance 
for the first waveguide mode of propagation [Black- 
band, 1963]. The displacement of the solid from 
the dotted line is a measure of the effect of the higher 
order waveguide modes which only propagate with 
appreciable amplitude to distances of up to about 
4500 km from the transmitter. The slope of the 
dotted line represents a change in velocity of about 1 
in 400 from day to night. 



3. Diurnal Change in Transmission Time 
Over Arctic Paths 

The details of the diurnal change in transmission 
time brought about by the change in ionospheric 
height are illustrated in figure 4 in which the upper 

1 This curve shows the change for the total sky wave. In practical applications 
the large values of At at short ranges are not apparent, being masked by the 
ground ray. 



diagram shows on a polar gnomonic projection the 
North Pole, and Arctic Circle. On midsummer day 
the line AB dividing the dark and light hemispheres 
of the earth's surface will be tangential to the Arctic 
Circle, and will appear to roll around this as the 
earth rotates. Consider a transmitter sited at T, 
and a receiver at R 1 . The inset (a) shows the 
diurnal variation observed at R u transmission time 
being plotted against local time at T. Starting at 
midnight with the whole path TR X in darkness and the 
apparent height of the ionosphere at its nighttime 
value, the delay is at its greatest. It will continue at 
this value as AB rolls around until the time D T 
when it is dawn at T. The effective height of the 
ionosphere then begins to decrease progressively 
along the path until D R when it is dawn at R lm 
This decrease is accompanied by a reduction in 
delay. Through the day the transmission time 
remains constant until sunset at R which occurs at 
a time S R , after this darkness spreads north until at 
S T the whole path is in darkness and the delay 
returns to its midnight value. 

This diurnal change exhibits the "diurnal change 
trapezium" discussed in section 2. However, more 
complex forms of pattern are possible. Some of 
these are shown in insets (b) , (c) , and (d) . 

Consider a receiver at R 2 sited so that dawn 
reaches both T and R 2 at the same time. As D T 
and D R occur at the same time the left-hand side of 
the pattern in the corresponding inset (b) shows a 
steep step, and the right-hand side is correspond- 
ingly reduced in slope. It is to be noted that at high 
latitudes the change from (a) to (b) can result from 
only a small change in the position of the receiver, 
such as from R x to R 2 . 
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Figure 3. Variation in diurnal change of transmission time with distance from the transmitter. 
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Insets (a) to) (c)fcO show the delay A t plotted against local time at T for 
the receiving sites Rj, R 2 , R 3 and R 4 . D and S are the times of dawn 
and sunset at T, D and S refer to the appropriate receiving sites. 
The lines NNand LLshow the steady values corresponding to the path all 
dark or all 1 i ght. 

Figure 4. Diurnal variations of delay in arctic regions. 

For a receiver at R 3 within the Arctic Circle there 
will be no time at which the whole path is in darkness 
and a curve of the form of inset (c) will result. It 
should be noted that the slow change in the fraction 
of path in darkness as AB rolls almost at right angles 
to the join of TR 3 , accounts for the rounded top to 
the pattern in the midnight hours. The time of 
maximum change is at the time when AB is tangential 
to the circle (i.e., midnight) at the point Q where 
the path TR 3 cuts the circle. The time of midnight 
at Q is calculable from its longitude. 

A more complex pattern is given for a receiver at 
i? 4 so sited that the path TR± crosses and recrosses 
the Arctic Circle. For such a path there are two 
periods each day when there is daylight for all the 
path. These are between dawn at T and sunset at 
R 4 and between dawn at i? 4 and sunset at T. Be- 
tween times there is darkness for part of the path 
and the delay increases. If the latitudes of T and 
R A are the same the delays in each period of partial 
darkness will be the same, otherwise that period 
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Figure 5. Phase patterns at various seaso?is. 



corresponding to darkness at the lower latitude end 
of the path will have the greater delay. 

The path between the transmitter at Seattle 
(NPG) and the receiving site at Farnborough 
crosses and recrosses the Arctic Circle in the same 
way as TR A of figure 4. 

The above explanation has been given in terms of 
ground sunrise and sunset largely because these are 
more definite than the more accurate concepts of the 
times when, looking from the daytime reflection 
height, the sun is seen to rise above or set below the 
upper boundary of the ozonosphere. Defining sun- 
rise and sunset in this way for reflection heights of 
67 km and ozonosphere extending to 35 km the mid- 
night sun at the reflecting layer would be seen as far 
south as 60°51' N. 

This effect could be allowed for in figure 3 by con- 
sidering the line AB as rolling round the circle of 60° 
51' N latitude instead of the Arctic Circle. It is to 
be noted that this latitude is that of the Shetland 
Isles, the southern tip of Greenland or Anchorage in 
Alaska and that for the purpose of this discussion of 
VLF propagation the "Arctic" is not confined within 
the Arctic Circle. 

4. Seasonal Variation in Diurnal 
Phase Change 

Consider an E-W path at middle latitudes as that 
discussed in section 2. With the change in the sea- 
sons D T and S T (and likewise D R and S R ) will Lome 
closer together for the short days of winter and 
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separate for the long summer days. The left-hand 
set of diagrams of figure 5 shows the corresponding 
lengthening of the trapezium in going from midwinter 
through the equinox to midsummer. It is to be 
noted that the basic form of the change pattern 
does not vary with season, and that At remains 
substantially constant. 

Figure 4 is drawn for midsummer day; as the year 
progresses the circle to which AB is tangential 
shrinks and, by the autumnal equinox, condenses on 
the pole. The line AB then crosses the pole, which 
is in darkness until the vernal equinox, and pro- 
gresses so as to reach the Arctic Circle from the other 
side at the winter solstice. Consider the effect of 
the shrinking of the circle to which AB is tangential. 
In the case of the path TR 3 as soon as the circle 
shrinks north of R 2 , the diurnal change of delay will 
be modified from the form of figure 4c to that of 
inset (a). Likewise as the circle shrinks so as not 
to cut the line TR±, there will be but one period of 
the day when there is daylight for the whole length 
of the path, and furthermore, there will be a period 
when there is darkness over the whole path— with 
these changes the diurnal change pattern for a path 
corresponding to figure 4d will become similar to 
that of figure 4a. When the line AB approaches the 
Arctic Circle at midwinter there will be a state w^hen 
there will be two periods each day when there will be 
darkness over the whole path and two when the 
path will be part dark part light, but no period of 
daylight for the whole path. This is complementary 
to the midsummer state in which there is no period 
for which the whole path is in darkness. 

The diurnal change patterns corresponding to 
such a path as TR± are shown in the right hand 
column of figure 5. It will be noticed that, unlike 
those of the left-hand column, these patterns change 
in form with season, and that At cannot have its full 
value during periods when the path is not both all 
dark and all light on the same day. The curves for 
this case illustrate the variation along the path 
Farnborough-Seattle. 

5. Measured Phase Change Patterns 

All the measurements recorded here have been 
made during the midsummer period when the 
peculiarities of the Arctic paths would be most 
marked. 

The paths monitored were 



(a) GBR (Rugby) 

(b) NPG (Seattle) 

(c) NPG (Seattle) 

(d) NPM (Hawaii) 



16.0 kc/s to 
Sondrestrom 

18.6 kc/s to 
Farnborough 

18.6 kc/s to 
Bod0 

19.8 kc/s to 
Bod0 



June 1961 
May 1961 
June 1961 
June 1961 



Station 


GBIl 


NPG 


NPM 


Bod0 


Sondre- 
strom 


Farn- 




borough 


Latitude 

Longitude 


52°22 / N 
1°11' W 


47° 3' N 
122° 20' W 


21°19' N 
157°50' W 


67°18' N 
14°26' E 


67°00' N 
50° 59' W 


51°17' N 
0°46' W 



The results of the measurements are summarized 
in the following table: 

Table 1 



Path 


Path 
length 


Length of 
of path 
having 

darkness 


Estimated 
At 


Measured 
At 


GBR — Sondrestrom 


3116 

7682 
7682 
6799 
10150 


1500 
2580 
2040 
1640 
4670 


12.0 
20.0 
15.8 
13.0 
36.4 


17 


NPG— Farnborough _. . 


25.9 


NPG— Farnborough 

NPG— Bod0 _ 


14* 
15.3 


NPM— Bod0 


34 



Distances expressed in km, At expressed in /xsec. 

In table 1 the second column shows the length of 
the path between transmitter and receiver. During 
the midsummer period, part of this path is in day- 
light at all times, and part will have times of light 
and dark. The length of the part changing from 
night to day on the day of measurement is shown in 
the third column. For the path NPG — Farnborough 
there are two entries the first corresponding to day- 
night changes at the western and the second corre- 
sponding to those at the eastern end of the path. The 
lengths quoted in column 3 are calculated on the 
basis of sunrise and sunset at an ionospheric height 
of 67 km, the sun being taken to rise or set not at the 
ground horizon, but at that of an ozonosphere 
extending to a height of 35 km. 2 

Of the five paths in the table only the first is short 
enough for the propagation of other than the first 
waveguide mode. For the other four the appro- 
priate values of At can be estimated by reference to 
the straight line of figure 3. These estimated values 
of At shown in column 4 of the table agree reasonably 
well with the measured values in the last column. 
The greatest difference is 5.9 /xsec of which some 
part is likely to be due to experimental error. 

For the shortest path, GBR-Sondrestrom, the path 
length was not great enough to justify the assump- 
tion that only the first mode was present with ap- 
preciable amplitude. However the estimated value 
of 12 //sec was based on this assumption — the dis- 
crepancy between this and the measured value of 
17 Msec probably arises from the presence of the 



2 The two figures of 67 and 35 km assumed for this calculation were based upon 
the following consideration. Measurerr ents [Bracewell, 1951], of daytime effec- 
tive height for VLF reflection made at near vertical incidence have given values 
of about 72 km. Because the effective height at glancing angle is expected to be 
somewhat less than this, and in any case the effective height is somewhat greater 
than the true height, 67 km seemed a reasonable estimate of the true height for 
reflection at large angles of incidence. The upper height of the ozonosphere is not 
definite but the value of 35 km was based upon the surveys of measured data 
published by Mitra [1952] and by D iitsch [1960] . The important parameter is the 
difference between these heights, that is the height of the reflecting layer above the 
ozonosphere. The value of 32 km which results from the above assumptions leads 
to reasonable agreement between theoretical and experimental values of At. 
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DEI AllS OF TRANSMISSION PATH SEATTLE - FAR MB OR OU OH 



SEATTLE 
LAT 4-7*03'N 
LONG \22'20\N 



FARNBOROUGH 
LAT 51*17' N 
LOWS 0*4b W 




Figure 6. Diurnal pattern for path Seattle — Farnboroagh. 
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Figure 7. Diurnal change curve for NPM at Bod<f>. 



higher modes. It is to be expected that their pres- 
ence would tend to increase At above the value ap- 
propriate to the propagation of a single mode. 

Some of the observed diurnal change curves are 
given in figures 6 through 8. For the path NPG- 
Farnborough a typical curve is shown in figure 6. 
This has two peaks in it, one corresponding to dark- 
ness at each end in turn. The curve has the general 
shape of that shown in figure 5 for Arctic summer. 
On the diagram are marked the times of high altitude 
sunrise and sunset (D T and S T ). Also marked are 
the calculated changes (as given in column 4 of the 



table) 
atP). 



and times of maximum change (midnight 



Similar data are given in figures 7 and 8 which 
relate respectively to the paths NPM-Bod0 and 
NPG-Bod0. 

6. Conclusions 

(a) The diurnal changes in transmission time for 

arctic paths differ from those expected on 
paths at middle latitudes. 

(b) That these changes vary in form and ampli- 

tude with season. 

(c) That these effects can extend to a latitude of 

about 60° 51' N. 

(d) That the seasonal variations can be estimated 

by simple methods. 
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Figure 8. Diurnal change curve for NPG at Bod<j>. 
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